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Abstract. We describe the spectra of neutral pions stemming from AuAu collisions at 
y/s = 200 AGeV in a ‘soft -I- hard’ model. The model is based on the assumptions that 
hadrons stemming from hard processes are described via perturbative quantum chromodynamics 
improved parton model calculations, while those stemming from the Quark-Gluon Plasma (we 
refer to as soft yields) can be described in a super-statistical model induced by multiplicity 
fluctuations. The obtained dependence of the parameters of the model on the event centrality 
is compared to what is observed in PbPb collisions at y/s = 2.76 ATeV. 


1. Introduction 

Among others, the phenomenon of jet-quenching supports that in ultra-relativistic heavy-ion 
collisions (AA) a new hadron source the Quark-Gluon Plasma (QGP) emerges besides the hard 
processes present in proton-proton (pp) and proton-nucleus (pA) collisions. We refer to hadrons 
stemming from these two types of sources as ‘soft’ and ‘hard’ yields. While perturbative 
quantum-chromodynamics improved parton model (pQCD) calculations [H El E] describe the 
spectra of hadrons stemming from hard processes (pp, pA collisions and the high transverse 
momentum pr ^ 6 GeV/c part of the spectra in AA collisions), pQGD methods fail to describe 
the Pt ^ 6 GeV/c part of hadron spectra in AA collisions. 

Disentangling soft and hard yields is a cumbersome problem for there may be low-px hadrons 
stemming from hard processes as well as moderately high-p^ ones created by soft processes |34j . 
Another difficulty is that due to the break-down of perturbation theory at low pT, we might 
not trust pQCD results at low momenta. In this paper, we demonstrate two possible ways of 
obtaining the low-pj’ part of the hard hadron yields: (i), we simply take the pQCD results; 
(ii), we use a function htted to the high-p^ part of the pQCD spectra to make an extrapolation 
to low pt- In the second case, we use a cut-power law function for the extrapolation, as such 
function turned out to provide an adequate description for hard hadron yields produced in pp 
collisions m-m- 

Once hard hadron yields (obtained in the (i) and (ii) scenarios) are subtracted from the 
measured transverse hadron spectra, we describe the remaining soft yields in the often used 
super-statistical framework [T2 [191 EOl EH EH EH EH SH jH]- this paper, we use the 
interpretation in which, the cut-power law shape of the soft hadron distribution originates from 
multiplicity fluctuations rather than from non-extensive entropy or transport [291E2] . 


The next section contains fits to neutral pion spectra measured in various centrality AuAu 
collisions at ^/s = 200 AGeV by the PHENIX Collaboration [l2]. To draw conclusions, we also 
discuss results obtained examining ALICE and CMS measurements [Odj . 

We note that similar soft and hard type separation of hadron yields has been suggested in 
models [221 [231 El] in which, spectra measured in pp collisions have been used as hard yields, and 
it has been conjectured that hard yields are suppressed at low px- Similar analysis’ of the soft 
part of hadron spectra without the subtraction of hard yields are also presented in the literature 

m-m- 


2. Transverse spectrum in the ‘soft + hard’ model 

As argued in Sec. [Hand in |34] . we make out the transverse spectrum of hadrons from the sum 
of hard and soft yields 
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As can be seen in the top-right panel of Fig. [H pQCD results describe transverse spectra of 
neutral pions stemming from various centrality AuAu collisions at y/s = 200 AGeV for pt > 6 
GeV/c. In the following, we examine two scenarios: 


i, We use pQCD results as hard yields for px < Q GeV/c too, and conjecture that the soft 
yields can be described by a cut-power law function. 
a, We fit pQCD results for px > Q GeV/c with one cut-power law. We use this fitted function 
to extrapolate the hard yields to low px- Finally, we use a cut-power law with different 
parameters to fit the soft part of the spectra. 


In this approach, both the px < Q GeV/c and the px > Q GeV/c regions contain soft as well as 
hard yields. The cut-power law finctions are 
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where i = soft or hard, m = m^o. The transverse mass is mx = y Px ^ isotropic 

part of the transverse flow with y* = 1/ — v‘f (see [M] for a detailed discussion). These yields 

have maxima at = 'jimvi. As long as these maxima are below the measurement range, 
which is the case in this analysis, the transverse flow, Vi cannot be determined accurately. 
As in this paper we analyse neutral pions, the arguments of /j-s may be approximated by 
[7j(mT — Vipx) — m]/Ti « px/T^’^ with the Doppler-shifted parameters 
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Results in the (z) scenario 

Fit results of Eqs. ([I|)-([2|) in the two scenarios explained above, are shown in Fig. [TJ The 
dependence of the obtained q and parameters of the soft and hard yields on the event 

centrality (number of participating nucleons Npart) can be fitted by (see Fig. Ej) 

Qi — Q2,i T part/‘2‘') , 

tDopp ^ T^ ^ + Ti\n{Npart) . (4) 

Though the actual value of the transverse flow velocity cannot be determined in this model 
from the spectra of charged hadrons, it may be guessed using the value of the QGP-hadronic 
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Figure 1. Top-left, transverse spectra of neutral pions stemming from various centrality AuAu 
collisions at y/s = 200 AGeV [l2]. Topmost graphs are central events, bottommost graphs are 
peripheral events in all panels. Solid curves; fits of Eqs. ©-(El), in which, both the soft and 
hard yields are htted cut-power law functions as described in the {i) scenario. Bottom-left, 
Ratio of measured data and fits in the (i) scenario. Top-right, ratio of measured yields and 
pQCD results. Bottom-right, ratio of measured yields and model results, in which, the hard 
yields are pQCD results and the soft yields are fitted cut-power law functions as described in 
the (ii) scenario. 




















































Figure 2. Fitted parameters of Eqs. © - (ED in the (i) scenario to PHENIX data on the 
spectra of vr'^-s stemming from AuAu collisions at y/s = 200 AGeV (see Eig. [1]). Eor comparison, 
results of fits to charged hadron spectra in PbPb collisions at y/s = 2.76 ATeV measured by the 
CMS and ALICE collaborations are also shown (see [S])- In both panels, data points belonging 
to soft hadron yields are colored in shades of blue, whereas those belonging to hard yields are 
reddish. Curves are fits of Eq. ©• 


matter transition temperature obtained from lattice-QCD calculations. As the values of fitted 
scatter around ~200 MeV at RHIC (and ~340 MeV at LHC), in case of a flow velocity 
of Vsoft ~ 0.5 at RHIC (and Vsoft ~ 0.6 at LHC), the real Tgoft values would scatter around 170 
MeV, which is close to the lattice result obtained e.g. in [43] . 

Qsoft and Tsoft may have statistical meaning. If in each nucleus-nucleus collision, hadrons 
stemming from soft processes have either Boltzmann-Gibbs or micro-canonical distribution, 
however, hadron multiplicity fluctuates according to the negative-binomial distribution, then 
the cut-power law function Eq. ([2]) emerges as the average hadron spectrum da [IB EE EB EEl 
EBlMlllol HI]. If 5 —)• 1, Eq. ([2]) tends to the Boltzmann distribution, while the multiplicity 
distribution tends to the Poissonian distribution. Such tendency can be seen in the left panel of 
Eig. E We may interpret this phenomenon as the soft yields are closer to the thermal, if either 
Apart or ^/s grows. As can be seen in the top panels of Eig. [T] and right panel of Eig. E the more 
central the collision, the bigger and hotter the soft yield. 

The behaviour of the parameters of the hard yields have to be determined by scale evolution 
equations of QCD in yet unexamined ways. Nevertheless, it can be seen that as y/s increases, 
Qhard increases as well, while the dimensionless parameter Thard/'/s appearing in the cut-power- 
law-type fragmentation functions introduced in [3], decreases. These trends seen in Eig. E are 
similar to what have been found in [3|- 















Q2,soft 

Q.2,,hard 

l^soft 

l^hard 

CMS 

1.058 ± 0.025 

1.136 ± 0.001 

-0.008 ± 0.005 

0.005 ± 0.0003 

ALICE 

1.074 ± 0.018 

1.131 ± 0.002 

-0.009 ± 0.004 

0.006 ± 0.0006 

PHENIX 

1.073 ± 0.016 

1.100 ± 0.002 

-0.005 ± 0.004 

0.000 ± 0.0006 


[MeV] 

rj^hard 

Tsoft [MeV] 

Thard [MeV] 

CMS 

310 ± 20 

126 ± 5 

9.9 ± 3.7 

5.3 ± 0.8 

ALICE 

266 ± 16 

194 ± 2 

11.5 ± 2.9 

-12.5 ± 0.5 

PHENIX 

165 ± 26 

192 ± 20 

9.3 ± 5.5 

18.7 ± 4.6 


Table 1. Parameters of Eq. Q obtained from fits shown in Fig. [2j 


Results in the {ii) scenario 

In the second scenario, the hard yields were taken directly from the pQCD calculations, and the 
soft yields were fitted by the cut-power law function Eq. ([2j) (see right panels of Fig. [1]). The 
obtained soft distributions were practically exponentials, thus we may conclude that if pQCD 
results mm are trustworthy at px < G GeV/c, the soft yields are likely to be thermal. This 
scenario has been thoroughly discussed in the literature [35]-[39]. 
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